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Abstract 

An experimental study of the normalized three-jet rate of b quark events with 
respect to light quarks events (light= £ = u,d, s) has been performed using the 
Cambridge and Durham jet algorithms. The data used were collected by the 
Delphi experiment at LEP on the Z peak from 1994 to 2000. The results are 
found to agree with theoretical predictions treating mass corrections at next- 
to-leading order. Measurements of the b quark mass have also been performed 
for both the b pole mass: and the b running mass: m^iMz)- Data are found 
to be better described when using the running mass. The measurement yields: 



m b {M z ) = 2.85 ± 0.18 (stat) ± 0.13 (exp) ± 0.19 (had) ± 0.12 (theo) GeV/c 2 
for the Cambridge algorithm. 

This result is the most precise measurement of the b mass derived from a high 
energy process. When compared to other b mass determinations by experi- 
ments at lower energy scales, this value agrees with the prediction of Quantum 
Chromodynamics for the energy evolution of the running mass. The mass 
measurement is equivalent to a test of the flavour independence of the strong 
coupling constant with an accuracy of 7°/ 00 . 
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1 Introduction 

In e + e~ collisions, data collected at high energies are predominantly of hadronic nature 
showing a multi-jet final state topology. At LEP, these data have led to measurements 
of many of the Standard Model (SM) parameters and allowed limits to be set on new 
physics processes. In some cases the original quark flavour and its mass have not been a 
critical issue in performing the measurement and therefore approximations using massless 
quarks have been sufficient for the required precision. This is well justified for inclusive- 
type observables like total cross-sections for which the correction due to massive quarks 
depends on the quark mass, m q , and on the energy of the process, Q, as m 2 c 4 /Q 2 . For 
b quarks (m q = m b ~ 3 — 5 GeV/c 2 ) at LEP I centre-of-mass energies (Q = M z c 2 ) this 
represents an effect of less than three per mille. On the other hand, for more exclusive 
observables, such as multi-jet cross-sections, the mass dependence transforms into terms 
proportional to m^c 4 / (Q 2 ■ y c ) where y c , the jet resolution variable, usually takes values 
much lower than unity, enabling these mass effects to become sizeable [1]. Differences in 
the multi-jet production rate for massive b quarks with respect to the corresponding rate 
for massless I quarks (£ = u, d, s) can then be as large as 3% to 20% for three- or four-jet 
final states. 

The large volume of data collected by the LEP experiments and the highly-effective 
techniques developed to identify the quark flavour of the jets, have increased the ex- 
perimental sensitivity for observables where the quark flavour is relevant [2]. Conse- 
quently theoretical predictions including mass corrections have become necessary to reach 
a proper understanding of such observables in order to interpret them as standard or new 
physics. In some cases these mass effects have only been computed at leading order (LO) 
but for some event-shape variables and in particular for the three-jet rate, calculations 
including next-to- leading order (NLO) terms also exist [3-6]. Rather precise experimental 
studies on the production of multi-jet events initiated by b quarks have then been allowed 
and performed at LEP and SLC [7-1 1] . Results obtained agree well with the predictions 
of the Standard Model, i.e., Quantum Chromodynamics (QCD). The quantification of 
these mass effects has allowed a verification of the flavour independence of the strong cou- 
pling constant to a precision of less than 1% and an extraction of a value of the b mass 
at the energy scale of the Z boson mass, Q = Mzc 2 , within an uncertainty of ~ 0.4 — 0.5 
GeV/c 2 . These measurements have in addition provided the first evidence of the running 
of the b mass, i.e. the evolution of this parameter as a function of the energy scale when 
compared to the values obtained from processes occuring at lower energy scales. 

A reduction of the uncertainty of the b mass determination could be accomplished by a 
combination of the present individual LEP and SLC results. This is certainly difficult and, 
at the end may improve the experimental precision only slightly, because the dominant 
errors are of systematic nature and common to all measurements. Hence, the best way 
to significantly increase the accuracy for this parameter is by a deeper understanding of 
the physics processes and correction procedures involved in the analyses. 

1.1 The b quark mass and the observable 

The b quark mass is a free parameter in the SM Lagrangian and therefore needs to 
be measured experimentally. Precise determinations of this parameter are very inter- 
esting both as a fundamental parameter and also to constrain models beyond the SM. 
Unfortunately, confinement of quarks inside the observed hadrons introduces additional 
complications not present in mass determinations of free particles as for instance leptons. 
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Quark masses need then to be defined within a theoretical convention and can only be 
inferred indirectly through their influence on hadronic observables. 

Among the different quark mass definitions, the most commonly used in high energy 
processes are the pole mass M b and the running mass m b {Q). The former is defined as 
the pole of the renormalized quark propagator and is gauge and scheme independent. 
The latter corresponds to the renormalized mass in the MS scheme and depends on 
the process energy (Q). These mass definitions are related to each other [12] and NLO 
calculations are needed in order to distinguish between the two. Physics is independent of 
the mass definition. However, when using perturbation theory at a fixed order to extract 
physics results a dependence appears on the mass definition as well as on the arbitrary 
renormalization scale \i. This is due to neglected perturbative and non-perturbative 
higher order terms and therefore it is possible than one scheme might be more convenient 
than another for a given purpose [13,14]. Specially because the relation between the 
two mass definitions does not have a good convergence behaviour due to renormalon 
ambiguities [15]. 

In the case of the running mass m b {Q), the largest part of its running occurs at 
low energy scales up to Mz/2. The exact mass running represents a basic constraint to 
theories beyond the SM, such as those implying the unification of the b quark and r lepton 
masses at the Grand Unified Theory (GUT) scale. The b quark mass has also important 
implications on Higgs searches since the partial decay width of the Higgs boson into b 
quarks is proportional to the b quark mass squared. In this case, it can be shown that 
the mass definition is also relevant for the accuracy of the theoretical prediction [16-18]. 

At low energy, the b quark mass is established from the measured spectra of hadronic 
bound states or the moments of the spectrum of the B decay products making use of 
non-perturbative techniques such as the Heavy Quark Effective Theory (HQET), Non- 
Relativistic QCD (NRQCD), QCD sum rules or lattice QCD. An attempt to average these 
determinations properly is presented in [13] obtaining the value m^m;,) = 4.24 ± 0.11 
GeV/c 2 . 

At high energy the b quark mass has been extracted from data collected at the Z 
peak at LEP and SLC. The first measurement of m b (M z ) was performed by the Delphi 
experiment with data collected from 1992 to 1994 at y/s ~ Mz- The observable used in 
this analysis was: 

3 {Vc) ~ Ri{y c ) ~ r§,-(vc)/i* (1) 

with r|^(y c ) and T q being, respectively, the three-jet and total decay widths of the Z into 
qq, where q = b or £ (£ = u, d, s quarks). The flavour q of the hadronic event was defined 
as that of the quarks coupling to the Z and the Durham algorithm was used for the jet 
clustering. The measured observable was compared with the theoretical computations 
of [3] and m b in the MS scheme was found to be [7] : 

m b (M z ) = 2.67 ± 0.25 (stat) ± 0.34 (had) ± 0.27 (theo) GeV/c 2 (2) 

where the quoted error was mainly due to the hadronization uncertainty. 

In this paper a new analysis to measure the b mass, performed with data collected by 
Delphi from 1994 to 2000 is presented. The data taken in the years before have not 
been considered because the Vertex Detector layout was changed in 1994, improving its 
capability since then (see Section 2). The same observable as used in the previous Delphi 
measurement, R^ e , has been used with two jet-clustering algorithms, Cambridge [19] 
and Durham [20]. The Cambridge algorithm has the advantage of having a smaller 
theoretical uncertainty [6]. A detailed study of how mass effects and the hadronization 
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process are implemented in the fragmentation models has led to a better control of the 
hadronization correction. The effect of the gluon-splitting rates into b and c quarks on 
the flavour tagging has also been taken into account. 

2 The Delphi detector 

A brief description of the most relevant components of the DELPHI detector for this 
analysis is given here. A detailed description of its design and performance can be found 
in [21,22]. 

Delphi was a hermetic detector with a superconducting solenoid providing a uniform 
magnetic field of 1.23 T parallel to the beam axis throughout the central tracking device 
volume. 

The tracking system consisted of a silicon Vertex Detector (VD), a jet chamber In- 
ner Detector (ID) and a Time Projection Chamber (TPC) which constituted the main 
tracking device in Delphi. At a larger distance from the interaction point the tracking 
was complemented by a drift chamber Outer Detector (OD) covering the barrel region 
(40° < 9 < 140°) and two sets of drift chambers, FCA and FCB, located in the endcaps. 

The VD was made of three coaxial cylindrical layers of silicon strips. From 1994 
onwards the outer and innermost layers were equiped with doubled-sided detectors with 
orthogonal strips, allowing the measurement of both Rip and z coordinates. In 1996, the 
VD was doubled in length and in 1997 a Very Forward Detector consisting of ministrips 
and pixels was added. Earlier DELPHI data taken in periods with a less complete VD 
setup are not included in this analysis whereas data collected, later, in the period of LEP2, 
1996-2000, which corresponded to the calibration runs at the centre-of-mass energy of the 
Z peak are used. 

Electron and photon identification was provided mainly by the electromagnetic 
calorimeter which was composed of a High Density Projection Chamber (HPC) installed 
inside the coil in the barrel region and a lead-glass calorimeter (FEMC) in the forward 
region. 

In order to measure the charged and neutral hadronic energy, Delphi also included the 
hadron calorimeter (HCAL), an iron-gas sampling detector incorporated in the magnet 
yoke. 

3 Hadronization correction of 

The hadronization correction to R^ e , i.e. the ratio of the observable at parton over 
hadron level, was computed in [7] using the string-fragmentation and cluster models 
implemented in Jetset 7.3 [23] and Herwig 5.8 [24] respectively, previously tuned to 
Delphi data [25] . Uncertainties coming from the tuned parameters of Jetset and from 
the different predictions of the two fragmentation models were taken into account, the 
latter being the highest contribution to the total error. 

In the present analysis more recent versions of the generators in which b quark mass 
effects are better modelled (Pythia 6.131 [26] and Herwig 6.1 [27]) were used. For 
the case of Pythia, different fragmentation functions were considered: Peterson [28] and 
Bowler [29] 1 . The final hadronization correction applied to our observable was the one 
given by Pythia 6.131 with the Peterson b fragmentation function since this model gives 

lr The Pythia parameters used are PARJ(55)=-0.00284 for Peterson and MSTJ(11)=5, PARJ(46)=1 and PARJ(47)=0.95 
for Bowler. 
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the best overall description of the Delphi data and the other two cases were only used 
to evaluate the model uncertainty. 

The model uncertainty was reduced to a negligible effect by performing the measure- 
ment in a restricted region of the phase space. Other sources of uncertainties such as the 
effect of the b quark mass parameter used in the generator were studied in detail and were 
shown also to be important. These two questions are discussed in the following sections. 

3.1 Hadronization model uncertainty 

The hadronization model uncertainty, a mo d, was evaluated as the standard deviation 
of the hadronization corrections predicted by the cluster model implemented in Herwig 
and the string-fragmentation model of PYTHIA using two heavy quark fragmentation 
functions, Peterson and Bowler. 

The hadronization corrections were found to depend on the S-hadron scaled energy, 
x| = 2E B _ haidTon /M z , the distribution of which is shown in Figure 1. As this quantity 
and the corresponding jet energy including the 5-hadron are highly correlated a new 
variable was defined instead: the 6-jet scaled energy x^(jet) = 2£ , 6_j et /Mz where -Eb_j e t is 
the energy of the jets originated by the primary b quarks in a Z — > bb event. The study of 
the dependence of the hadronization corrections with x^(jet) led to the conclusion that 
if the cut x^(jet) > 0.55 is applied to both b quark jets the model uncertainty is reduced 
by a factor of 4 (see the right hand plot in Figure 1 and Figure 2 ) . 




Figure 1: x| distribution for the different Monte Carlo generators (left). Hadronization 
model uncertainty (i.e. standard deviation of the hadronization corrections predicted 
by Herwig and Pythia with Peterson and Bowler heavy fragmentation functions) as 
a function of the mean of the x^(jet) that is defined in the text and for y c = 0.01 for 
Cambridge (right). 
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Figure 2: Hadronization model uncertainty as a function of the y c for the CAMBRIDGE 
(left) and DURHAM (right) algorithms. 



3.2 The b mass parameter in the generator 

In this section, the effect of the b quark mass parameter used in Pythia on the 
hadronization correction is discussed. The result of this study also applies to other 
generators which contain similar features. 

In order to describe b quark mass effects, Pythia uses a set of three b quark mass 
parameters: the kinematical mass, Mff me , used in the parton shower (PS) process, the 
constituent mass, M£ onst , used during the hadronization process and finally the known 
B hadron masses. The constituent mass is also used to derive masses for predicted but 
not yet observed B hadrons. In the model these three masses are not connected to each 
other and, as a consequence, mass effects at parton level do not automatically propagate 
to the hadronization process, as they physically should. This results in a dependence of 
the hadronization correction on the b quark mass of the generator. 

If the various mass parameters are connected by, for instance, making the constituent 
and kinematical b quark masses equal to each other and by deriving all B hadron masses 
from the corresponding quark masses using the hadron mass formula [30], this depen- 
dence of the hadronization correction factor is completely removed. This feature of the 
generators was also noticed in previous studies [7-1 1] even though the exact cause of this 
behaviour was not identified. Unfortunately, this argument cannot define the value to be 
assigned to the mass parameters of the standard generator in which the quark masses are 
not connected. For that purpose, the following procedure was applied: 

• In order to assess the precision of the massive calculation implicit in the parton 
shower generator, its prediction for R h ^ and that of the NLO calculation was first 
compared. The method was to change the input mass in the NLO calculations 
to minimize their overall difference. Then the difference of the input mass values 
was evaluated. For the parton shower the so called kinematical mass, M£ tne , was 
employed and for the NLO calculations the two mass definitions were considered. In 
the case of the running mass the corresponding value was transformed to the pole 
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mass 2 . The value A(M fe fcme — M b ) ~ 15 MeV/c 2 was obtained in the case of the 
running mass and A(M^ me — M b ) ~ 500 MeV/c 2 for the pole mass. These differences 
were later considered as the uncertainty associated to the effective mass definition 
of the parton shower. 

• The values of the b quark mass measurements determined from low energy processes 
were then used as input to the mass parameter, M^ me , of the generator entering 
in the hadronization process. In order to select the mass value to be used in the 
present analysis, various possibilites were explored. A direct determination of M b 
from reference [31] gave M b = 4.98 ± 0.13 GeV/c 2 . A second possibility is to use 
the average value for the running mass calculated in [13] as: m b {m b ) = 4.24 ± 0.11 
GeV/c 2 , which could be transformed into a pole mass value of M b = 4.99 ± 0.13 
GeV/c 2 . A third value is also available using Delphi data from the semileptonic 
B decays for which the relevant scale is that of the B hadron masses, leading to 
M b = 5.00 ± 0.16 GeV/c 2 [34]. All these results are compatible and have similar 
accuracy. The value which was used in the generator to compute the hadronization 
correction was that obtained as the average of all low energy measurements [13]: 

m b (m b ) = 4.24 ±0.11 GeV/c 2 , or, M b = 4.99 ± 0.13 GeV/c 2 . (3) 

For cross-check purposes, values of the b quark mass were also extracted using Delphi 
data alone with the modified generator for which the set of the three b quark masses are 
connected to each other. Two different observables were employed for this study: the 2/32 
distribution 3 of b over £ events normalized to the total number of b and £ events, R be (y32) 
and the minimum angle between b quark and gluon jets when every event is forced to 
three jets. Both quantities are correlated with the observable i?| used to measure the b 
quark mass and therefore their role in the present analysis is limited to qualitative checks. 
The first observable gave a fitted value for the b quark mass of the modified generator of 
M b = 4.93 ±0.13 GeV/c 2 and the second one M b = 4.95 ±0.11 GeV/c 2 . These results are 
thus consistent with the choice of the mass parameter and the above quoted uncertainty. 

4 Experimental determination of Rf 

First the sample of Z hadronic decays, i.e. Z — > qq events was selected. Then the b 
and £ quark-initiated events were separated using the Delphi flavour tagging methods 
and later a cut on the b quark jet energy was also performed in order to discard those 
events with large hadronization correction (see Section 3.1). 

The jet-clustering algorithms Cambridge and Durham were applied to both tagged 
samples to obtain the R}§ observable at detector level. Data were then corrected for 
detector and tagging effects and for the hadronization process to bring the observable to 
parton level. 

4.1 Event selection 

The selection of Z hadronic events was done in three steps (as in [7]): 

• particle selection: Charged and neutral particles were selected in order to ensure a 
reliable determination of their momenta and energies by applying the cuts listed in 
Table 1; 

2 The 3-loop relationship between m h (m b ) and M b [32] with a s (M z ) = 0.1183 ± 0.0027 [33] was used. 
3 j/32 is the y c transition value in which a 3-jet event becomes a 2-jet 
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• event selection: Z —> qq events were selected by imposing the global event conditions 
of Table 1; 

• kinematic selection: In order to reduce particle losses and imperfect energy- 
momentum assignment to jets in the selected hadronic events, further kinematical 
cuts were applied. Each event was clustered into three jets by the jet-clustering algo- 
rithm (Cambridge and Durham) using all selected charged and neutral particles. 
The cuts of Table 1 were then applied. 

After applying these cuts to the data a sample of 1.4 x 10 6 (1.3 x 10 6 ) hadronic Z 
decays was selected for the Cambridge (Durham) algorithm. 



Charged 

Particle 

Selection 


p > 0.1 GeV/c 
25° < 9 < 155° 
L> 50 cm 

d < 5 cm in R<p plane 
d < 10 cm in z direction 


Neutral 
Cluster 
Selection 


E > 0.5 GeV, 40° < 9 < 140° (HPC) 

E > 0.5 GeV, 8°(144°) < 9 < 36°(172°) (FEMC) 

E > 1 GeV, 10° < 9 < 170° (HCAL) 


Event 
Selection 


N ch >5 
E ch > 15 GeV 
YliQil < 6, i = 1, ...,N ch 
No charged particle with p > 40 GeV/c 
45° < 9 thrust < 135° 


Kinematic 
Selection 


Nf l > 1 per jet 

Ej > 1 GeV, j = 1,2,3 

25° < 9j < 155°, j = 1,2,3 

Planarity cut: J^a faj ^ 359 °' i < h hj = 1 ? 2 , 3 



Table 1: Particle and hadronic-event selection cuts; p is the particle momentum, 9 the 
particle polar angle and 9thrust the thrust polar angle (with respect to the beam axis in 
both cases), L the measured track length, d the closest distance to the interaction point, 
qi the particle charge, E the cluster energy, N c h the number of charged particles, and 
E c h the total charged-particle energy in the event. The kinematic selection is based on 
the properties of the events when clustered into three jets by the jet algorithm. Nj h is 
the jet charged multiplicity, Ej the jet energy, 9j the jet polar angle and 0^- the angular 
separation between the pair of jets ij. 



4.2 Flavour tagging 

The b and light (£ = u, d, s) quark-initiated events need to be identified. Delphi 
has developed two different algorithms for b tagging based on those properties of B 
hadrons that differ from those of other particles: the impact parameter [35] and the 
combined technique [36]. The former makes use of the most important property for the 
selection of B hadrons, their long lifetime, and discriminates the flavour of the event by 
calculating the probability, P£, of having all particles compatible with being generated 
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at the event interaction point. The second technique, besides the impact parameter 
of charged particles, uses other discriminating variables: the transverse momentum of 
any identified energetic lepton with respect to the jet direction and, in case a secondary 
vertex is found, the total invariant mass, the fraction of energy, the transverse momentum 
and the rapidities of the charged tracks belonging to the secondary vertex. An optimal 
combination of this set of variables defined for each reconstructed jet is performed, leading 
to a single variable per event, X e ff ev . When the aim was to measured i?g for Cambridge, 
this was also the algorithm used to compute the combined tagging variable, and the same 
for the case of Durham. 

Figure 3 shows the distributions of Pg and X e ff ev obtained for the selected real and 
simulated sample of Z hadronic decays. For the case of the simulated data, the contri- 
bution of each quark flavour is also indicated. 




Figure 3: Event distribution of P£ (left) and X e ff ev (right) when Durham is used to 
form jets. The real (points) and simulated (histogram) data are compared. The specific 
contribution of each quark flavour is displayed as derived from the Delphi simulated 
data. The cuts used to tag the b and I quark sample are also indicated. 

Taking into account the stability of the final result (see Figure 4 left), the impact 
parameter method was used for I tagging by imposing P£ > 0.07. The resulting purity 
of the sample and efficiency of the selection were Pg = 82% and q = 51%, respectively. 
For b tagging both techniques were observed to be equally stable. The combined method 
was used requiring X e ff ev > —0.15 since higher purities could be reached for the same 
efficiency The final purity of the sample and the total efficiency were = 86% and 
tb = 47%, respectively, where this efficiency value also takes into account the hadronic 
selection. 

In order to perform the cut on the b quark energies (see Section 3.1), an identification 
of the gluon and b quark jets was required for 6-tagged events. The two tagging techniques 
can also provide a discriminant variable per jet and therefore both are available for jet 
identification. Again, based on a stability argument (see Figure 4 right), the combined 
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Figure 4: (Left) R% = r| (y c ) /T q at parton level as a function of the purity of the g-tagged 
sample, P q , for q — b,£, when Cambridge is used to form jets. (Right) -Rjj as a function 
of the jet identification purity for the Durham algorithm. 



technique was used to identify the pair of jets most likely to come from b quarks by 
requiring X e jf ev > —0.5 (where now X e ff ev is only computed with the tracks contained 
in the pair of jets which gives the maximum X e ff ev ). This results in a b jet purity of 81% 
per jet in each event and a tag efficiency of 90%. 

Once the b quark jets were identified their energy was computed from the jet directions 
using energy-momentum conservation and assuming massless kinematics. Figure 5 shows 
the x|.(jet) distribution for real and simulated data. The cut Xg(jet) > 0.55 was then 
applied for both b jets. The purity and contamination factors of the b and ^-tagged 
samples obtained after the cut are shown in Table 2. 



Method 


Type q 


£ -> q (%) 


c -> q (%) 


b -> q (%) 


Impact Parameter 


i 


82 


15 


3 


Combined 


b 


2 


7 


91 



Table 2: Flavour composition of the samples tagged as £ or b quark events. (£, c, b) — > q 
refers to the fraction of true q' events in the g-tagged sample. 



4.3 Data correction 

Once the b and £ quark hadronic events were selected from the collected data, the 
jet-clustering algorithm was applied to get the Rl§ observable at detector level, R^ e ~ det . 
In order to bring the observable to parton level the method of the previous Delphi 
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Figure 5: Xg(jet) distribution for real and simulated data for three-jet ft-tagged events at 
y c = 0.015 for Cambridge. 

measurement was used [7]. A correction to obtain pure b and £-quark samples was applied 
in this procedure and the flavour composition uncertainties were accounted for by the 
tagging uncertainty. 

The Delphi full simulation (Delsim), which uses Jetset 7.3 to generate the events 
that go through the detector simulation, was used to compute the detector correction. 
A reweighting of the events was done in order to reproduce the measured heavy quark 
gluon-splitting rates [37] (g cZ = 0.0296 ± 0.0038 and g bl = 0.00254 ± 0.00051) in the 
simulation. 

A recent version of Pythia 6.131, tuned to Delphi data [25] and with the kinematical 
b quark mass parameter set to M b = 4.99±0.13 GeV/c 2 , was used to get the hadronization 
correction. 

The magnitude of the detector and hadronization corrections for i?| are shown in 
Figure 6. At the y c value chosen for the final result (y c = 0.0085 and y c = 0.02 for Cam- 
bridge and Durham, respectively) the detector correction is about -2.5% for Durham 
and 5% for CAMBRIDGE. The hadronization correction is 1% for CAMBRIDGE and half 
as big for Durham. 

4.4 Experimental uncertainties 

Apart from the statistical uncertainties, different sources of systematic uncertainties 
were considered. They can be divided into two groups: those due to the hadronization 
correction and those due to the detector correction. 

• hadronization: 

The following sources of uncertainties in the hadronization correction have been 
taken into account: 
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Figure 6: Detector and hadronization corrections applied to the measured R£ for Cam- 
bridge and Durham. The detector correction, C det , brings the observable to hadron 
level, and the hadronization correction, C had , brings it from this stage to parton level. 

— uncertainty in the tuned parameters of Pythia that are relevant in the frag- 
mentation process. This contribution was evaluated by varying these parameters 
(Aqcd, <5o; e b, a ) within ±1 standard deviation around their tuned central 
values, taking into account correlations [25]; 

— uncertainty due to the choice of the hadronization model. It was calculated as 
the standard deviation of the three hadronization models used (see Section 3.1); 

— uncertainty from varying the value of the b quark mass parameter in the gen- 
erator within the error of 0.13 GeV/c 2 about its chosen central value of = 
4.99 GeV/c 2 (see Section 3.2). 

• detector: 

The uncertainties in the detector correction, including selection efficiencies, accep- 
tance effects and the tagging procedure, are due to imperfections in the physics and 
detector modelling provided in the simulation. The following sources were consid- 
ered: 

— gluon-splitting: The measured cc and bb gluon-splitting rates were varied within 
their uncertainty and the effect on the measured observable was taken as the 
gluon-splitting error; 

— tagging: The related uncertainty was evaluated by varying the tagging and 
mis-tagging efficiencies within their uncertainties: Ae^/e^ = 3% and Ae e b /e e b = 
Ae£/e£ = 8% evaluated as in [38] (being e q q the fraction of q' tagged events in 
the true g-quark sample). For this purpose, £ tagging was considered equivalent 
to anti-6 tagging, i.e. Ae 9 e = Ael for q = b,c,£ for the same cut value; 

— jet identification: The cut applied to distinguish the b quark jets from the gluon 
jet in a b tagged event was varied in order to obtain cut efficiencies (i.e. the 
fraction of events which pass the cut applied to select the two b quark jets among 
the 3 jets) between 80% and almost 100%. Half of the full variation observed 
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in the measured R% at parton level was taken as the uncertainty due to the jet 
identification. 

4.5 Results for R\ at hadron and parton level 

Figure 7 shows, as a function of the y c , the measured i?jf corrected to hadron level to- 
gether with the curves predicted by the Pythia and Herwig generators. The statistical- 
only and total uncertainties can also be seen in this figure. The Pythia curves are shown 
independently for the Peterson and Bowler b fragmentation functions. For large values of 
y c , both generators describe the data well. The measured R^ e ratio and its uncertainties 
are also presented in Tables 3 and 4. 




Figure 7: i?| at hadron level as a function of y c compared with Pythia 6.131 (with 
Peterson and Bowler fragmentation b functions) and Herwig 6.1 predictions, using the 
Cambridge (left) and Durham (right) jet-clustering algorithms. 

The result for R^ e obtained at parton level is shown in Figure 8 as a function of 
y c together with its statistical and total uncertainties. The LO and NLO theoretical 
predictions in terms of the pole and running masses (M& = 4.99 GeV/c 2 and rrib(Mz) = 
2.91 GeV/c 2 ) are also shown in the plot. In the case of Cambridge the LO prediction is 
already reproducing the measured data, indicating a better convergence in the theoretical 
calculations than for Durham. The results for the individual years of data taking are 
compatible (see Figure 9). 



5 Comparison with NLO massive calculations 

The measurement of the R% observable at parton level obtained in the previous section, 
when compared with the NLO massive calculations of [3,6], can be used either to extract 
the b quark mass assuming ot s universality or to test a s flavour independence taking the 
b quark mass measured at threshold as an input. 
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Figure 8: R% as a function of y c obtained at parton level compared with the LO and 
NLO theoretical predictions calculated in terms of a pole mass of Mb = 4.99 GeV/c 2 and 
in terms of a running mass of m^iMz) — 2.91 GeV/c 2 . 
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Figure 9: E!§ at parton level obtained for each analysed year for a fixed y c for Cambridge 
(left) and Durham (right). The error bars represent the statistical error. The vertical 
lines show the average value with its statistical and total error. The y 2 per degree of 
freedom of the average is 0.7 and 1.2 for Cambridge and Durham respectively. 
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Vc 


T)bt—had 
R 3 


&stat 


® 9~ splitting 


O'tag 


& jet— id 


0.005 


0.9617 


±0.0034 


±0.0003 


±0.0014 


±0.0016 


0.01 


0.9544 


±0.0044 


±0.0010 


±0.0025 


±0.0019 


0.015 


0.9560 


±0.0052 


±0.0014 


±0.0031 


±0.0021 


0.02 


0.9632 


±0.0059 


±0.0018 


±0.0036 


±0.0024 


0.025 


0.9639 


±0.0067 


±0.0021 


±0.0039 


±0.0026 


0.03 


0.9629 


±0.0074 


±0.0024 


±0.0041 


±0.0029 



Table 3: Rf at hadron level at different y c with jets reconstructed with Cambridge. 



Ve 


j^bi—had 


®stat 


® g— splitting 


O'tag 


& jet-id 


0.005 


1.0194 


±0.0033 


±0.0002 


±0.0011 


±0.0008 


0.01 


0.9690 


±0.0039 


±0.0007 


±0.0025 


±0.0010 


0.015 


0.9613 


±0.0047 


±0.0012 


±0.0031 


±0.0012 


0.02 


0.9583 


±0.0056 


±0.0016 


±0.0036 


±0.0014 


0.025 


0.9596 


±0.0062 


±0.0018 


±0.0040 


±0.0015 


0.03 


0.9611 


±0.0070 


±0.0022 


±0.0043 


±0.0017 


0.035 


0.9606 


±0.0076 


±0.0025 


±0.0045 


±0.0019 


0.04 


0.9630 


±0.0083 


±0.0027 


±0.0046 


±0.0021 


0.045 


0.9626 


±0.0089 


±0.0029 


±0.0048 


±0.0022 


0.05 


0.9687 


±0.0097 


±0.0032 


±0.0050 


±0.0024 



Table 4: Rf at hadron level at different y c with jets reconstructed with Durham. 

5.1 Determination of the b quark mass 

In order to extract the b quark mass from the experimentally measured a value 
of y c must be chosen for both Cambridge and Durham jet algorithms. The value used 
was that which gave the smallest overall uncertainty on the measurement while staying 
in the region where the hadronization correction remains flat. In this way it was also 
guaranteed to keep far enough from the four-jet region. The selected values found to best 
fulfill these requirements were y c = 0.0085 and y c = 0.02 for Cambridge and Durham, 
respectively, where the four-jet rates are 4-5% and 2-3% in each case. 

The b quark pole mass, M , could be extracted from the measured R^ e using the NLO 
expression of i?| in terms of M& [3,6]. Theoretical sources of uncertainty were the \x 
scale dependence, the identification of the b quark mass parameter in the generator (see 
Section 3.2) and a s . 

The measured b quark pole mass was found to be, 

M b = 4.19 ± 0.23 (stat) ± 0.17 (exp) ± 0.25 (had)+°;™ (theo) GeV/c 2 (4) 
when Cambridge is used to reconstruct jets with y c = 0.0085 and, 

M b = 4.47 ± 0.31 (stat) ± 0.24 (exp) ± 0.24 (had)^ (theo) GeV/c 2 (5) 

when Durham is used instead with y c = 0.02. Although compatible within errors, these 
values are low compared with the results obtained when the b pole mass is measured at low 
energy (as for example 4.98 ± 0.13 GeV/c 2 [31]). The measurement error is dominated 
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by the uncertainty from the identification of the b quark mass parameter in the generator 
with the b pole mass which contributes to the theoretical error. 

The running mass was also obtained using the NLO computations of i?|f from refer- 
ences [3,6], in this case, in terms of the running mass at the M z energy scale: m b (M z ). 
The theoretical uncertainty was estimated by considering the following sources: 

• dependence on the renormalization scale: The /i scale in the theoretical expressions 
was varied from Mz/2 to 2Mz and half of the difference between the result obtained 
on rrib(Mz) was taken as the /i scale error; 

• mass ambiguity: Starting from the NLO calculation of R^ e in terms of the pole mass 
M b , the value of M b could be extracted and transformed to m b (M b ) which was later 
evolved to m b (M z ) by means of the Renormalization Group Equations. This is also 
a valid procedure to extract m b (Mz). At infinite orders the result derived in this 
way and that obtained directly from the original NLO calculation in terms of the 
running mass should be the same. The difference between the results obtained from 
the two procedures was then also considered as a conservative indication of the size 
of the unknown higher order corrections; 

• a s : a s (M z ) = 0.1183 ± 0.0027 [33] was varied within its uncertainty. The spread of 
values obtained for m b [Mz) was considered as the error due to the a s uncertainty. 

The results obtained for m b {Mz) were, 

m b (M z ) = 2.96 ± 0.18 (stat) ± 0.13 (exp) ± 0.19 (had)l^ (theo) GeV/c 2 (6) 

when Cambridge was used to reconstruct jets and, 

m b [M z ) = 3.42 ± 0.25 (stat) ± 0.18 (exp) ± 0.20 (had)l^ (theo) GeV/c 2 (7) 

in the case Durham was the algorithm employed. 

The theoretical uncertainty expressed in this way is highly asymmetric due to the 
mass ambiguity. Hence the interval covered by the extreme values of the theoretical 
uncertainty originating from this mass ambiguity was considered as the whole range of 
theoretical uncertainty and the measurement of m b (M z ) was set to the mean value of 
this region. The effect on the mass value is a shift in the order of ~ -100 (-200) MeV/c 2 
for Cambridge (Durham). The same criteria were also adopted in previous work [7,9] 
and in the present case leads to: 

m b [Mz) = 2.85 ± 0.18 (stat) ± 0.13 (exp) ± 0.19 (had) ± 0.12 (theo) GeV/c 2 (8) 

when Cambridge was used to reconstruct jets and, 

m b [Mz) = 3.20 ± 0.25 (stat) ± 0.18 (exp) ± 0.20 (had) ± 0.24 (theo) GeV/c 2 (9) 

if the Durham algorithm was used. 

The contribution of the individual uncertainties is given in Table 5. The result obtained 
with Cambridge is more precise than the one obtained with Durham mainly because 
of the smaller theoretical uncertainty, leading to a total error of ±0.32 GeV/c 2 instead of 
±0.44 GeV/c 2 . 

5.2 Test of a s flavour independence 

The measurement of R^f can alternatively be used to test a s flavour independence 
exploiting the relation introduced in [7]: 
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aj/aj = R be - H(m b (M z )) + A • - H(m b (M z )) - 1), (10) 

71 

where H(rrib(Mz)) is the theoretical mass correction and the factor A depends on the jet 
reconstruction algorithm and y c , taking values between 2 and 6 for all possible circum- 
stances of the present analysis. 

Taking the average b quark mass from low energy measurements, m b {m b ) = 4. 24 ±0.11 
GeV/c 2 [13], as the input b mass value, the ratio a b s /a e s is found to be, 

a b J a e s = 0.999 ± 0.004 (stat) ± 0.005 (syst) ± 0.003 (theo) (11) 
for Cambridge and 

a b J a e s = 0.990 ± 0.006 (stat) ± 0.006 (syst) ± 0.005 (theo) (12) 
for Durham. These results verify a s universality at a precision level of 7-9°/ 00 . 



6 Conclusions and discussion 

A new determination of the b quark mass at the M z scale has been performed with the 
Delphi detector at LEP. The same observable as for the previous Delphi measurement 
[7] was studied, now also using the Cambridge jet clustering algorithm in addition to 
Durham. The results obtained with Cambridge for m b (M z ) were found to be more 
precise, giving: 

m b (M z ) = 2.85 ± 0.32 GeV/c 2 . (13) 

This constitutes a substantial improvement with respect to the previous Delphi mea- 
surement in which m b (M z ) was determined to be 2.67 ±0.50 GeV/c 2 . This is mainly due 
to the improved evaluation of systematic errors as has been described in this paper. 

When using the theoretical prediction of Rf for the Cambridge algorithm the data 
are reasonably well described by the theoretical calculation, already at leading order, 
using the value m b (Mz) = 2.91 GeV/c 2 inferred from the low energy measurements (see 
Figure 8). The higher-order terms contributing to the calculation of the observable appear 
to already be accounted for in the running of the mass and therefore a faster convergence 
seems to be achieved in comparison with the b pole mass. However for DURHAM the 
situation and the behaviour are different as in fact both theoretical predictions at LO 
are equally distant from the data using both mass definitions and NLO calculations are 
certainly needed to describe the data. The value for the b pole mass determined in this 
case was: 

M b = 4.1918;™ GeV/c 2 . (14) 

The present measurement has been performed in a restricted region of the phase space 
to have a better control of the fragmentation process. New versions of the generators, 
Pythia 6.131 and Herwig 6.1, where mass effects are much better reproduced, have 
been used to correct the data. 

The study of the way mass effects are implemented in the generators, described in 
Section 3.2, has led to a more reliable hadronization correction. The pole mass definition 
was shown to be the one to be used in the generator and the uncertainty of this identi- 
fication on the present analysis has been quantified. It constitutes the dominant source 
of the present error. The effect of the bb and cc gluon-splitting rate uncertainties of the 
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Monte Carlo on the detector correction has also been taken into account. The observable 
B!§ is also presented at hadron level for different y c values in view of future versions of 
the generators with a better understanding of the hadronization process which could then 
allow for an improved measurement. 

The result obtained by this analysis with Cambridge for m b {Mz) is shown in Figure 
10, together with other LEP and SLC determinations at the M z scale. It is compatible 
with the other measurements and is the most precise. The data collected by Delphi 
have also been used to determine the b quark mass at a lower energy scale near threshold 
using semileptonic B decays [34]. The value obtained in that analysis is also shown. 
The difference between the two measurements is significantly larger than the overall 
uncertainty: 

A(m b (m b ) - m b (M z )) = 1.41 ± 0.36 GeV/c 2 . (15) 

Hence, for the first time, the same experimental data allow values for the b quark mass 
to be extracted at two different energy scales. The results obtained agree with the QCD 
expectation when using the Renormalization Group Equation predictions at the two 
relevant energy scales of the processes involved. These observations together with the 
average value of the b quark mass determinations at threshold [13], m^m^), are shown 
at their corresponding scales in Figure 10. 

Alternatively, universality of the strong coupling constant has also been verified with 
a precision of 7°/ 00 . 

For data combination purposes, the above results supersede the previous DELPHI 
measurements on this subject [7]. 
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Figure 10: The evolution of m^Q) as a function of the energy scale Q. The m^Mz) 
measured by LEP and SLC are displayed together with their total and statistical errors. 
The shaded area corresponds to the band associated to m^Q) when evolving the average 
value obtained at m^rrib) [13] up to the Mz scale using the QCD Renormalization Group 
Equations with a s (Mz) = 0.1183 ± 0.0027. All these measurements are performed at the 
Mz energy scale but for display reasons they are plotted at different scales. The result 
obtained using Delphi data at low energy from semileptonic B decays [34] is also shown. 
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Cambridge 


j^b£—had 


■pybt—part 


m b (M z ) 




(y c = 0.0085) 


(y c = 0.0085) 


GeV/c 2 


Value 


0.9527 


0.9646 


2.85 


Statistical Data 


±0.0033 


±0.0034 


±0.14 


Statistical Simulation 


±0.0024 


±0.0025 


±0.11 


Total statistical 


±0.0041 


±0.0042 


±0.18 


Fragmentation Tuning 





±0.0010 


±0.04 


Fragmentation Model 





±0.0025 


±0.11 


Mass parameter 


— 


±0.0036 


±0.16 


Total hadronization 




±0.0045 


±0.19 


Gluon-Splitting 


±0.0008 


±0.0008 


±0.03 


Tagging 


±0.0022 


±0.0021 


±0.09 


Jet identification 


±0.0018 


±0.0020 


±0.09 


Total experimental 


±0.0030 


±0.0030 


±0.13 


Mass Ambiguity 






±0.11 


//-scale (0.5 < fi/M z < 2) 


— 


— 


±0.04 


a s (M z ) 


— 


— 


±0.01 


Total theoretical 


— 


— 


±0.12 


Durham 


T)bi—had 


-p>bt—part 


m b (M z ) 




(y c = 0.02) 


(y c = 0.02) 


GeV/c 2 


Value 


0.9583 


0.9626 


3.20 


Statistical Data 


±0.0045 


±0.0045 


±0.20 


Statistical Simulation 


±0.0033 


±0.0033 


±0.15 


Total statistical 


±0.0056 


±0.0056 


±0.25 


Fragmentation Tuning 





±0.0015 


±0.07 


Fragmentation Model 





±0.0022 


±0.10 


Mass parameter 


— 


±0.0034 


±0.15 


Total hadronization 




±0.0042 


±0.20 


Gluon-Splitting 


±0.0016 


±0.0016 


±0.07 


Tagging 


±0.0036 


±0.0035 


±0.15 


Jet identification 


±0.0014 


±0.0018 


±0.08 


Total experimental 


±0.0041 


±0.0042 


±0.18 


Mass Ambiguity 






±0.22 


//-scale (0.5 < fi/M z < 2) 






±0.10 


a s {M z ) 






±0.02 


Total theoretical 






±0.24 



Table 5: Values of at hadron and parton level and of m^M^) obtained with Cam- 
bridge and Durham algorithms and the break-down of their associated errors (statis- 
tical and systematic) for y c = 0.0085 and y c = 0.02 respectively. 
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